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Abstract 


Quantum walks are roughly analogous to classical random walks, and like classical walks 
they have been used to find new (quantum) algorithms. When studying the behavior of 
large graphs or combinations of graphs it is useful to find the response of a subgraph to 
signals of different frequencies. In so doing we can replace an entire subgraph with a single 
vertex with frequency dependent scattering coefficients. 

In this paper a simple technique for quickly finding the scattering coefficients of any 
quantum graph will be presented. These scattering coefficients can be expressed entirely in 
terms of the characteristic polynomial of the graph’s time step operator. Moreover, with 
these in hand we can easily derive the “impulse response” which is the key to predicting 
the response of a graph to any signal. This gives us a powerful set of tools for rapidly 
understanding the behavior of graphs or for reducing a large graph into its constituent 
subgraphs regardless of how they are connected. 
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1 Introduction 


In classical computer science random walks have proven to be a useful tool for understand¬ 
ing and developing new algorithms and techniques. The same has been true for quantum 
walks [T] and quantum computers. A number of these quantum walks have already been 
experimentally implemented, some using trapped ions mm and others using photons in 
optical networks 03-m The goal of this paper is to provide a set of powerful and computa¬ 
tionally cheap tools for rapidly understanding the behavior of graphs in terms of scattering 
signals. 

A graph is a set of vertices with edges connecting those vertices. In a classical random 
walk we imagine a particle that inhabits a vertex and has some probability of moving to 
connected vertices. In a discrete random walk time is divided into integer steps and the 
probability of the particle jumping to an adjacent vertex in a given step is described by a 
stochastic matrix. Like classical walks, the time can be either continuous [8] or advance 
in discrete steps mm ■ But unlike classical walks, with quantum walks the time step 
operator is unitary rather than stochastic. 

In a classical walk the Hilbert space is composed entirely of the set of vertices, but 
in a quantum walk that isn’t sufficient. In a nutshell, there isn’t enough information in 
a vertex state alone for time-reversibility (an essential characteristic of unitary processes) 
because at the very least the particle needs to “know” where it was in the previous time 
step. In this paper we’ll use “edge states” to more elegantly encode this information m- 
For example, the edge state | A, B) is the state on the edge between vertices A and B that 
points from A to B. This is exactly equivalent to a particle on vertex B that was previously 
on vertex A. In this edge state formalism each vertex hosts a unitary operator that takes 
all of the incoming states and maps them to outgoing states. 

Those already familiar with quantum walks will probably be more familiar with the 
“coin space” [lj formalism, which is more directly analogous to classical random walks; the 
Hilbert space of discrete time coined quantum walks is the tensor product of the position 
space and the coin space m- That is, each vertex is amended with an ancillary “coin 
space”. The coin keeps track of where the particle previously was, and for this reason the 
dimension of the coin space associated with a vertex is always greater than or equal to the 
degree of that vertex. 

Finally and most importantly, in a quantum walk the particle is in a superposition of 
states described by probability amplitudes (as opposed to probabilities). When a position 
measurement is made the probability of finding a particle in the state | ip) on the edge |e) 
is | (e|"0) | 2 , in adherence to Born’s rule. 

Several previous studies have investigated scattering on quantum graphs. Similar to 
the formality found in this paper, semi-infinite lines (“runways” of attached edges) are 
attached to a graph. On these runways the time step operator passes the particle from 
one edge to the next sequentially, either toward or away from a particular vertex in the 
graph, and in this way the particle enters and exits the graph. In [8j scattering theory was 
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used to show that tree graphs could be used to solve some kinds of decision problems using 
continuous time walks. A discrete time scattering theory approach was fleshed out in m 
where the connection between the connection between the number of steps to get through 
a graph and the transmission amplitude was found, as well as some more general results 
on the reflection and transmission amplitudes for a graph. In m it was shown how the 
scattering matrices for subgraphs could be used to construct the scattering matrix for the 
overall graph. 

In this paper we show that the semi-infinite runways of edges can be replaced with a 
single edge. This allows us to analyze problems using strictly finite graphs, removing the 
issue of non-normalizable states. More importantly, we show that the scattering coefficients, 
as well as the response to any incoming states, can be described entirely and succinctly 
by the characteristic polynomial of that finite graph’s time step operator. The time step 
operator is dictated by the structure of the graph, so this allows us to immediately see the 
relationship between the scattering coefficients of a graph and the structure of that graph. 

Specifically, we find that the scattering coefficients of a graph are determined by the 
eigenstates and eigenvalues of the time step operator (which is now a finite matrix). Fortu¬ 
nately, the techniques described in this paper do not at any time require these quantities to 
be calculated; instead we find that the characteristic polynomial itself is all that is needed. 

In section 2 we consider the case of a single runway. The problem and the exact 
definition of the effective reflection coefficient are defined rigorously. Having only a single 
runway leads to some surprisingly compact results that are explored here. We then use the 
frequency-dependent effective reflection coefficient to derive the graph’s impulse response. 
With this in hand we are able to rapidly calculate (with a convolution) the response to any 
arbitrary signal. 

In section 3 we address the challenges of using a reflected signal to gain information 
about a graph and a theorem is proven which describes this difficulty explicitly. Here we 
learn the relationship between the eigenvalues of the time step operator and the length 
of a signal on the runway necessary to detect the effect of those eigenvalues. This is an 
important tool for understanding the computational time of algorithms. 

In section 4 we explore the case of a graph attached to multiple runways, and the 
effective scattering coefficient between them. A very powerful theorem is derived that 
allows for the rapid and simultaneous calculation of each of these coefficients in terms of 
the resolvent. 

Finally, in the appendix (the second half of this paper) there are a series of examples 
that put all of the theorems and techniques described in this paper to use, demonstrating 
how simple they are in practice. 
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2 Basic Framework 


The situation in question is an arbitrary graph, G, attached to an infinite “runway” of 
edges. The vertices on the runway are labeled 0, 1, 2, ... where 0 is a given vertex of G. 
We define the unitary time step operator on the runway as the one that passively moves 
each edge state one step. I.e., U| j,j + 1} = |j + 1 ,j + 2) and U |j + 1, j) = \j, j — 1). The 
behavior of U in the graph G is not specified here. 

Our goal is to replace the graph G with a single vertex and to encode all of its behavior 
into that one vertex. A set of constant reflection/scattering coefficients doesn’t contain 
nearly enough information to simulate the behavior of a complicated graph, but if we allow 
them to be frequency-dependent, then this goal is attainable. 

The frequency-dependent scattering coefficient is defined such that it replaces the entire 
graph with a single reflection coefficient at vertex zero with value .S'(A). When dealing with 
multiple runways attached to the same graph we can consider them one at a time, since 
the time step operator is linear. 


X 2 X 

-« c > t < G 

2 \ 1 S(X) 1 S(\) 0 

Figure 1: A signal that advances by A every time step produces a “reflection” that is shifted 
by some phase. The scattering coefficient is defined such that the graph can be replaced 
with a single vertex that reflects with S'(A). 

In this section we’ll first consider a graph with one connected runway. It will be 
demonstrated that we can understand the response of a graph connected to an infinite 
runway by looking at the characteristic polynomial of the graph alone. 

We define an operator, U, on G and the runway such that U|l, 0) = | in) and U| out) = 
10,1), and Vn > 0, U|n — 1, n) = \n, n + 1) and U|n + 1, n) = |n, n — 1). 

An incoming pure momentum state takes the form ^ +l \j + 1 ■ j)■ If vertex 0 is 

completely reflective with reflection coefficient r, then the response is r ^ n ~' > \jij + 1) 

(see figure 1). Clearly, 


|T> = £ A J+1 \j + l,j) + + 1) (1) 

3=0 j =o 

is an eigenstate with eigenvalue A. 

In this form it’s easier to see how this is a signal and reflection. After n time steps 
this will take the form A n |T) = F^o^ n+J+1 |-? + kj) + J2'jLo r ^ n ~ : ’\jG + 1)> an d after 
n + 1 time steps the state is A n+1 |^) = A n+J+2 |j + 1 ,j) + Y^=Q r ^ n+1 ~^\ji j + !)■ 
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The coefficient of 10,1) is r times the coefficient of 11,0) in the previous time step, which 
is exactly as it should be. 

If instead of a single simply-reflecting vertex there is a graph attached to vertex 0, then 
the eigenstate now takes the form 

OO OO 

l*> = E Xj+1 \j + ^i) + E S ^ X ~ j \ 3,j + 1) + \G) (2) 

3=0 j =0 

where |G) is the component of the eigenstate contained in G. In this way we can 
define an “effective reflection coefficient”, 5(A). Unless G is a single vertex, 5(A) will be 
a non-constant function of A. In either case, equation [2] is a A-eigenstate of the graph and 
runway. 



ctC 


Figure 2: U and U a . By selecting the correct value of a we can produce a A-eigenstate 
of U a that is identical to the A-eigenstate of U on the edges they have in common (all of 
which are in G). 

In order to find 5(A) we create a new operator, U Q , that reflects back into the graph 
rather than transmitting into or receiving from the runway. That is: 

XJ a | out) = a | in) (3) 

We will find that there is a simple relationship between a and 5(A), and that we can 
determine the correct value of a by tuning it such that U a has A as an eigenvalue. 

First we need to derive a few properties of the characteristic polynomial of U Q . 
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Theorem 2.1. C(z, a) = |U Q — zl\ = b(z)(f(z) + ag(z)), where f(z), g(z), and b{z) are 
polynomials in z. f(z) and g(z) share no common roots, the roots of b(z) sit on the unit 
circle, and the roots of f(z) sit strictly within the unit circle. 

Proof 

This is easy to immediately verify by inspection of the matrix U a — zl. a appears once, 
so every term in the characteristic polynomial either contains an a or doesn’t. Clearly, the 
characteristic polynomial is affine in a. 

We can collect the terms with and without a’s into two polynomials. Trivially, those 
polynomials can be labeled b(z)f(z) and ab(z)g(z), where b(z) is the collection of all of 
the factors common to both polynomials. 

Since the roots of b(z) are independent of a , and since U Q can be unitary (when 
|a| = 1), the roots of b(z) are eigenvalues of a unitary matrix and therefore have modulus 

1. 

Clearly, f(z) = |Uo — zl\, where Uo := U Q | a=Q . Define |'I'o} = a\out) + |G) to be a 
normalized eigenstate of Uo and rj to be a root of f(z). Since Uo| out) = 0| in), we know 
that (m|\ko) = 0. When |a| = 1 we know that U a is unitary and therefore 
\h? = (vkoluJUol^o) 

= (d'o| ^Ul — a*\out)(in\J (U a — a\in)(out\) jd'o) 

= (\E'o|UaU a |\I/o) — 2 Re [a*('ko|ortt)(in|U Q ,|'I , o)] + |a| 2 o\out)(in\in)(out\'& o) 

= (d'ol^o) — 2 Re [a*a*(m|U a |^o)] + M 2 M 2 
= 1 — 2 Re [a*a* (aa)] + |a| 2 |a| 2 
= 1 — |a| 2 |a| 2 
= 1 — |a| 2 
< 1 

If a = 0, then |^o) is a bound eigenstate, and rj would actually be a root of b(z). 
Therefore all of the roots of f(z) are inside the unit circle. 

□ 

In everything that follows b(z) either doesn’t play a roll, isn’t relevant, or factors out. 
So, it will be suppressed. So far we’ve replaced one unknown variable, 5(A), with another, 
a, however this is a step forward because we can quickly find a closed solution for a. 

Theorem 2.2. 5(A) = \ = 

Proof 

This is the essential trick of this paper. 



The A-eigenstate ofU takes the form (T) = Yl'kLo A fc+1 |A: +1, k) + 5(A) J2T=o ^~ k \k, k + 
1) + A5(A)|oitf) + | in) + |G). 

The A-eigenstate of U a takes the form |'k a ) = 4| out) + | in) + |G). 

Restricted to G, these two states are the same. But whereas the coefficient of | in) in 
|\k) is dictated by the incoming signal on the runway, in | T Q ) it’s dictated by the coefficient 
of | out) and the value of a. By tuning a to the correct value we’re “feeding the output to 
the input” so that (m|U n |'L) = (m|U”|\I / a } = A n , Vro. 

Equating the coefficients of the | out) states in these two eigenstates we see immediately 
that 5(A) = 4. We can then solve for a using the characteristic equation. When A is 
an eigenvalue we have that 0 = /(A) + ag( A) and therefore a = — It follows that 

5(A) = — • Note that S(z ) is a meromorphic function of z. 

□ 


To reiterate and make clear what 5(A) is: so long as the runway and G are in the 
A-eigenstate, we can replace the G with a reflection coefficient, 5(A), at vertex 0. In this 
way we can describe the graph’s reaction to an infinite signal with frequency A = e ld . 


2.1 Particulars for a Single Runway 

In the case of a single input and output XJ a is unitary for |a| = 1. This unitarity has a lot 
of consequences, but in particular S(z ) takes the following form: 


S(z) = 1 
a 


9(z) 

m 


where \rjj\ < 1, V). 


ffo tt 1 Z7 lj 
2 s 11 2: — 77 7 

i J 


(4) 


In this section we’ll explore the special case of a single runway. 


Theorem 2.3. C(z,a ) = goaz d C* Q, 4) 5 Vz,Va / 0 where C* indicates the coefficients 
are conjugated, d is the degree ofC(z,a), and go is the constant term of g(z). Equivalently, 

f(z) = z s HU ( z ~ ’ d{z) = goz d f* (i) = go UU ( X “ zr] i) ’ where d ' + s = d - 

Proof 

In what follows assume that |a| = 1. This means that U a is unitary, and the roots of 
the associated characteristic polynomial, C(z, a), all have modulus 1. While this proof will 
only consider f(z) and g(z), it works in exactly the same way for b{z)f(z) and b(z)g(z). 

Let C(z,a) = IIfc=i ( z ~ h) = Efe=o fk zk + a Efc=o 9kZ k - Note that f d = 1 and 
9d = 9d- 1 = 0, since when the determinant was taken any term with a necessarily did not 
include 2 diagonal elements (2 powers of z). 
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/(Afc) + a5(A fc ) = 0 
<7(Afc, a) = 0 
^ (C(Afc, a))* = 0 


C*{X* k ,a*) = 0 


« C *(sd) =° 

« 0 = /• (A) + Is* 



Therefore, (7(2,a) and C*(-,-) 
az d C*(l^) = az d f* Q) + z d g* Q) 
same set of zeros. Therefore, az d C* 
bC(z, a) 


have the same set of zeros. It also follows that 
is a polynomial in 2 and a which, again, has the 
, and C(z,a) are proportional to each other. 


= ^ d c'(bi) 

= az d Et=o4*? + nSl;=o9fc? 

= *EL 0 r k z d - k + Et= 0 9iz d ~ k 

= a Xb=0 fd-k z>C + Xb=0 9d-k zk 


f bfk = g* d _ k 
1 bg k = f*_ k 
Now, 

fd = 1 =*> b = 

We now have that C(z, a) = goaz d C* Q, ^). 

Since the constant term in a characteristic equation is equal to the determinant, 1 = 
|a<7o| = bo| 5 which means that 0 is not a root of g(z). Because g ( i = gd-i = 0 we have that 
/o = fi = 0, which implies that s > 2 

Keep in mind that the result above is merely a statement about the polynomial C(z, a). 
It is true regardless of the value of a. 


(7(2, a) = g 0 az d C* (-, 

\z aJ 

in general, Va / 0. Or equivalently, 


(5) 


f(z) = g 0 z d g* 
g(z) = 90 z d r Q) 

f(z) and g(z) are said to be “reciprocal polynomials” of each other. 


( 6 ) 

(7) 
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□ 


We can say even more about the characteristic polynomial. The zeros of f(z) and g{z) 
have a very particular behavior and relationship. 

Theorem 2.4. C(z,a ) = T (z - rjj) +a g 0 n (l-zrij) 

3 3 

V ^ V V V 

/( 2 ) ff( 2 ) 

where 0 < \r)j\ < 1, Vj. 

Proof 

We already know that \r]j\ < 1 from theorem 
the last theorem. It follows that for iy / 0, f(rjj ) = 0 44 g j = 0. Therefore, if 

f(z) = z s n, (z - Vj), then g(z) oc ■ ( z ~ Fly (l ~ zr/* ) • With g 0 the constant 

term in g(z), we can write g(z) = Ft y (l — zr fj ) • 

Notice that f(z) and g(z) share a root if and only if rjj = 4r or 1 77 ^ | = 1. But this is 

rjj 

exactly what we expect for the roots of b{z). □ 

The above statement about the roots of f(z) applies more generally; that is, it continues 
to apply when there are multiple inputs and outputs. 


2.1 


and that f(z) = goz d g* (*) from 


Theorem 2.5. When |a| = 1, the solutions ofC{z,a ) are distinct. 

Proof 

When |a| = 1 we know that U a is unitary. An immediate consequence of which is the 
fact that U Q is diagonalizable and expressible as U Q = \ where Pa is a projection 

operator onto the A-eigenspace. Each of these projections can be expressed in terms of the 
resolvent, which in turn can be written as a power series in a — a 0 near ao, where ao is 
any arbitrary point on the unit circle. 

This implies that the projection operators can likewise be expressed as a power series in 
a — ao, and since Pa = | V\)(Pa 11 it follows that the eigenvectors share the same property. 
Finally, since U Q and |V\) are power series in a —ao, and U q |Va) = A|V\), we can see that 
the eigenvalues themselves, A, are power series in a — ao- 

Now define co(z)(z — Ao Y = f(z) + aog(z). Note that according to the last theorem 
g{z) 7 ^ 0 when \z\ = 1 . 
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o = /(A) + ag{\) 

= /(A) + a 0 g{\) + (a - a 0 )g( A) 

= c 0 (A)(A - A 0 ) 4 + (a - a 0 )g(X) 

=> (A - A 0 ) f = (a - «o) 

=r- A = Ao + O (-^cr — ao) 

However, since the eigenvalues are expressible as a power series in a — qq, t = 1. 
Therefore, because «o is arbitrary, the degree of any zero of f(z ) + ag(z ) is one when 
|a| = 1. 

□ 

In this proof it was important that |a| = 1 because it ensures that U a is unitary. For 
a finite set of values of a (off of the unit circle) we find that f(z) + ag(z ) can have higher 
degree roots, however at those points we find that U a is no longer diagonalizable and the 
degenerate eigenvalues correspond to generalized eigenvectors. 


Theorem 2.6. When a loops once around the unit circle the eigenvalues cyclicly permute 
one step. That is, looping a changes A j —>• Xj+i and A^y —> Ai, where arg (Ai) < arg (A 2 ) < 
• < arg (X d ). 


We know that looping a once (returning it to its original value) can’t change the 
spectrum of the eigenvalues, so the effect must be a permutation. In addition, since the 
eigenvalues are always distinct for every value of |a| = 1, this permutation must be cyclic 
(the eigenvalues can’t ”slide past each other” on the unit circle). 

So we know that looping a produces a permutation of the eigenvalues of the form 
A j —> X j+t (where A = Ao). The only question that remains is the value of t. 

Define A = e ie . The eigenvalues satisfy 
0 = /(A) + ag( A) 


=> ~f (e W ) = ag (e zd ) 

=> —e is6 n-Li {e ie - r,j) = ago Uf =1 (l - e%*) 

=> ~e lse n -Li {e 16 ~ Vj) = ag 0 Y[f=i e l ° (e l ° - iy)* 
=> -e' se uU ( el6 - *&) = <*90**° TlU ( eld ~ *Y 


=> lo g(«) = i7T + l°g(ffo) + i(s ~ d’)6 + Y?- =1 log 

=> i arg(a) = in + i arg(g 0 ) + i(s - d')6 + ^=1 ^ ar g ( e * 0 “ Vj) 
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We now have a relation between the zeros of f(z ) and the phase of a. 


d' 

arg(a) = n + &rg(g 0 ) + (s - d!)0 + 2 ^ arg (V 0 - r y^j ( 8 ) 

3 = 1 

At this point we allow 9 to smoothly increase by 2ir, then take the difference. Since 
\rjj\ < 1 , the angle between e* e and ly sweeps from 0 to 27 t monotonically. 

A arg ( a ) = (s — d')2ir + 2 i ^ 

A arg (a) = (s + (1')2tt 

A arg (a) = 2nd 

Looping a given eigenvalue once around the unit circle causes a to loop s + d' = d times. 
Looping an eigenvalue once is a permutation of the form Xj —>• \j+d = A j. It follows that 
if looping a once produces a permutation of the form A j —> \j+t, then looping A j means 
that a loops j times. But we know that looping an eigenvalue once requires a to loop d 
times, and therefore t = 1 . 

□ 


Theorem 2.7. Any eigenvalue A, such that |A| = 1, can be induced by choosing the correct 
value of a. Moreover, this value is unique. 

In the last theorem it was shown that the eigenvalues, which are functions of a, cyclicly 
permute when a loops around the unit circle. These functions are continuous, so every 
value between these eigenvalues exist for some value of a as well. Moreover, arg( A) is a 
strictly monotonic function of arg(a), and therefore the value of a is unique for a given 
eigenvalue. 

This can be seen by first showing that -ggarg (e ie — r/j > \ when \r/\ < 1. This can be 
proven by either using the inscribed angle theorem to establish a lower bound or by direct 
calculation. It follows that 

arg (a) = i r + (s - d')6 + i log (g 0 ) + 2 J2f=i arg (e ld - 7y) 

=> m ar g («) = {s-d!) + 2 J2jL 1 §$arg (e ld - ry) > (s - d!) + 2d' Q) = s > 0 
=> §garg (a) > 0 

Therefore arg(a) and arg(A) are strictly monotonic functions of each other. This mono¬ 
tonicity ensures the uniqueness of a for a given value of A by making sure that arg(a(A)) 
doesn’t double back on itself. 

□ 
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2.2 Arbitrary Inputs for a Single Runway 

In the language of signal analysis, the last section was a derivation of the “frequency 
response” of the graph G. We can define the input x[n] (output y[n}) as the amplitude on 
the state j 1,0) (state 10,1)) at time step n (time step n + 1). The input can be encoded 
onto the runway in an initial state of the form YlkLo + 1> k). 

At time step n, the overall state of the graph and runway will be: 

OO OO 

y~]s[fc + n\\k + 1 ,k) + '^2’y[n- k - l\\k + l,k) + |G) (9) 

k =0 k =0 

That is, we can describe the coefficients of the edge states on the runway with the “signal 
function” x[n] and the “response function” y[n] which are series of complex numbers, one 
for each integer n. 

At the nth time step x[n] is the coefficient of 11,0) and y[n — 1] is the coefficient of 
10,1). This is so defined such that if S(z) = r, where r is constant, then y[n\ = rx[n]. 
When S(z) is constant this holds true for any x[n], but for non-trivial S(z) we have 


x[n\ = z n ,\/n 


y[n] = S(z)x\n], Vn 


( 10 ) 


This is nothing more than a restatement of the definition of S(z), as described in eq. 
[2j Using S(z) we can derive an expression for the “impulse response”, h[n], which is 

( 1, n = 0 

the response, y\n], produced from the input x\n] = 5\n], where 5\n] = < ’ . _ is the 

[ 0, n 0 

Kronecker delta function. From the definitions of x[n], y[n ], and the impulse response itself 
we can say that 


h[n] = (0, l|U n+1 |l, 0) (11) 

That is; the impulse is just the state j 1,0) at time zero and the impulse response is just 
the coefficient of 10,1) read at each sequential time step. 

We already know that for a simple reflection (such that U|1,0) = r|0,1)) y[n\ = rx[n] 
and therefore h[n] = rd[n]. This can be seen from either of equations [lO] or [TlJ The response 
of any signal can be found using the fact that y[n] = (x*h)[n] = h[k\x[n — k\. So rather 
than being a single example, the impulse response is the key to finding the response to any 
input (15) . 


Theorem 2.8. If h[n\ is the impulse response (that is, y[n} = h[n] when x[n] = S[n]), then 
for any signal function x[n], y[n\ = (x * h)[n}. 

Proof 
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A graph’s response to signal functions is a map, T, from the space of sequences of 
complex numbers to itself. When we say that y[n\ is the response to x[n], we can write 
this more succinctly as T (x[n|) = y[n\. T inherits linearity from the linearity of U. 

Define x\~ = x[k\, \/k. We do this to distinguish between the function x[n] and the value 
of the function evaluated at k. 

We can use the Kronecker delta function to break x[n] apart into a sum of simple 

Xfc ,n = k 


signals. Trivially, for any fixed value of k, XkS[n — k] = XkS[k — n] = 
Therefore, if we sum over k we reconstruct the full function: 


0 ,n / k 


x[n] = ^2 XkS[n — k] 
k 

From this and the fact that T (<5[n]) = h[n] it follows that 
y[n\ 

= T ( x[n ]) 

= T('£ k x k 5[n- k ]) 

= J2k x kT'(S[n- k]) 

= J2k x kh[n - k] 

= (x * h ) [n] 

□ 


Theorem 2.9. The impulse response, h[n], is given by h[n] = 2 si|„| =1 ^ n 1 S(z)dz = 

_3_ X n-l 

2ni Z f(z) az 


Proof 

For a single frequency, x\n ] = A n , we find that 
V M = Efc h[k]x[n - k] 

= E k m \ n ~ k 
= E k h[k] X~ k X n 
= {Ek h i k \ x ~ k ) x [ n ] 


For any fixed value of A we already have a way of writing this (thm. 2.2). Therefore 


Rewriting A = e lu \ 


J>[fc]A- fc = S(A) 

k 


y( x ) 

/(A) 


( 12 ) 
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S(A) = £ fc M^]A- fc 

=> Efc h[k\e~ iuk = S ( e i0J ) 

=► J2k h[k]e iuJ ( n ~V = e inuj S ( e iu} ) 

=> / 0 2w Efc /#]e Mn_fc) ^ e inw 5 (e iw ) dw 

27r/i[n] = J Q 27r e ,nw 5 (e^) dw 
d[n] = ^ / 0 2?r e m0J S (e* w ) du 

d[n] = 2 ^ ,^|=i z n_1 S l (z) dz (z = e* w , dz = ie luJ du) 

=* %] = -2h i^=i zn ~ l m dz 

This is a perfectly nice integral, since the zeros of f(z ) are all on the interior of the 
unit disk. 

□ 


Theorem 2.10. For a single input and output, h[n} = flod[n — s] + J2j ^jVj 

~9oU k (^) ,3=0 

i-M 2 n ( Izmul \ ,-_^ n 

90 v s+1 IIMj V 'h-m ) ’d ^ U 


where Flj = < 


Proof 

From the previous theorem the impulse response is: 


1 


h[n] = - — 


2iri 


N=i 


z n ~ 1 ^f\dz = 

f(z) 2ttz 


z n ~ s ~i TT i—^ dz 
l*M \ z ~ Vk 


(13) 


Using residue calculus we can solve this directly 
h[n\ 


90_ £ z n- 

im !r\z\=a z 

— 'n* 

1_ ^ dz 
z-m 



don fc (4-) 

1 5[n — s 

\-aoJ2ji 1 - \vj\ 2 )Ylk^j 

V Uj-Vk ) 

i vr- 1 

don fc (4^) 

1 5[n — s 

l + u[ m l p n ¥l ( 

' 1 ~ r u r ii\ 

, Uj-Vk ) 

v n 


So the impulse response is a set of exponentially decaying signals corresponding to the 
zeros of /(z). 
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3 Quantum Sounding 


Since there is a closed form for the response to any signal that is a function only of the 
eigenvalues of Uo, we can (at most) find the spectrum of a graph’s eigenvalues. We are now 
equipped to ask the question “What can be learned about a graph attached to a runway 
by means of a signal on that runway?” and even ask the more practical question “How 
difficult is it to do so?”. 

The challenge we face is that the closer r/ is the unit circle, the more difficult it is to 
detect. From eq. |4j we know that (for a single runway) S(z) = — ~ S Elfc ~z^t ' 

= 1. 


l- 


z % 


z-Vk 


When \z\ = 1 we can easily verify that |S'(z)| = 1. More specifically, for each k, 

We’ll look at each of these individually and since each has modulus 1, we can concern 
ourselves entirely with their phase. 


Here we assume that 0 -C \r]k\ < 1. When z r/k we can see that = ~Vk 


z-r)k 


'%■ 


Importantly, 1 k is approximately constant for most values of z on the unit circle. 

z Vk 

rjki and therefore 


z-r/k 


This is because -V = 

Vk \Vk\ 

has one pole inside of the unit circle and one zero outside of it. As a result, when 

Z V k 

we apply the argument principle we find that if z runs around the unit circle in a positively 
oriented loop, then A arg ( 1 ^ fc ) = —27r. 

V z Vk J 


It follows then that a zero of f(z) has very little impact on S(z ) except when z is 
within a small neighborhood of that zero and within that small neighborhood the phase 
suddenly jumps by —27r. We’ll now make this a little more rigorous to find the extent of 
this neighborhood. 


Theorem 3.1. If ij = (1 — 5)e lT is a root of f(z) and x[n] = e md , we find that can only 
be detected when \9 — r| = 0(5). Moreover, the phase of the reflection coefficient decreases 
by 2 tt in this neighborhood of i 7 . 


Proof 


Without loss of generality, we can assume that rj = 1 — 5. We define = 

anc ^ quickly find that 


e* = 


-1 + 


<5 2 (1 + cos ( 9 )) 


+ i 


(—25 + 5 2 ) sin (9) 

2(1 — <5)(1 — cos (9)) + 5 2 


2(1 — <5)(1 — cos (9)) + 5 2 

Clearly, for 9 0, lim, 5 ^.o = —1- Define c = j? = 0(1). The imaginary part of this 

last equation for small values of 5 and 9 is 
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• / /\ (-25 + 5 2 ) sin (c<5) _ —2cc5 2 + 0(<5 3 ) _ 2c 

Sm[<p) ~ 2(1 — <5)(1 — cos (c<5)) + <5 2 c 2 (5 2 + <5 2 + 0(<5 3 ) “ “^Tl + W 

So the window for which \<j)\ < ^ is approximately —5 < 6 < 6 and the window for 
which \<f> — 7r| > d (for which is different from -1) is approximately — < 8 < ^<5- 

This is the statement of the theorem. 

□ 


In appendix |A.4 there is an example of how a graph structure can be used to detect 
the phase of a reflection coefficient. 


From theorem 2.1 we know that b(z)f(z) = |Uo — zl\, which means that the degree of 
the polynomial f(z) is less than or equal to the number of edge states in G. S(z) = — 
as established in theorem 2.2 Since the zeros of f(z ) are all within the unit circle and the 
zeros of g(z ) are all outside, we can apply the argument principle to S(z) around the unit 
circle to conclude that 


| A arg (£( 2 ))) < 2tt\G\ 

The equality is achieved if there are no bound states. So, with very little effort, we 
have an algorithm that provides a lower bound for the dimension of a graph’s Hilbert space 
(the number of edge states). 
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4 Scattering: Multiple Inputs and Outputs 

Define Uo to be the time step operator of a finite graph, G, with some states prepared as 
either loose inputs or outputs. An input state has no pre-image and an output state has 
no image. That is: Uo|out) = 0 and Uq l \in) = 0. 

If we wish to “splice” a runway onto the graph we first choose one input state, | in 3 ), and 
one output state, | out}.), and define as Utf*)|l,0> = | i U j) and U^| out k ) = |0,1). 

The behavior on the runway is the same regardless of which states on G it connects with, 
so Vj, k and n > 0, n — 1, n) = |n, n + 1) and U^|n + 1, n) = |n, n — 1). 

A A-eigenstate of U^ fc ) necessarily takes the form: 


|HO = ^2 |^A n+1 |n + l,n) + Sjk( A)A n \n,n+ 1} 

n=0 


+ | irij) + XSjk{X)\outk) H- 


(14) 


As before we introduce another operator, U^ = Uo + a\irij)(outk\, that reflects back 
into the graph rather than communicating with the runway. For both of these operators 
U^' fe) | outi) = UW| outi) = 0, Vi ^ k. 

(jJA 

A A-eigenstate of Uq ' necessarily takes the form: 

|T a ) = \irij) + ^\outk) 4- (15) 

If the A-eigenstates are identical on G, then clearly ^ = XSjk(X) and therefore Sjk( A) = 

a * 


Theorem 4.1. Sjk{ A) = — 9 fff , where Cjk(z) = 
characteristic polynomial of U« ^ . 

Proof If A is an eigenvalue of U a k \ then 0 = f(X)+agjk{X). By matching the coefficient 

of the \outk) state in |T a ) and IT) we find that Sjk( A) = ^ = — 9 j k (X) • 

□ 


U, 


ti k ) 


-z I 


= f(z) + agjk(z) is the 


This is essentially the same as thm. |2.2| with one unimportant difference. In this case 

f n (jk)\ 

is no longer unitary since Uq \outf) = 0 for £ k. As a result we can no longer say 
that |a| = 1, however this has no impact on the proof. In fact, it is to be expected that 

we see that for \z\ = 1 , Ylk \Sjk( z )\ 2 = 1 so 


4.3 


for multiple runways |a| > 1. From thm. 
we can conclude that | Sjk(z)\ < 1 and |a| = 


|5jfc(z)| 


> 1. 
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For any operator M we call (M — zl) 1 the “resolvent” of M. The resolvent has many 
fascinating properties [16] . and here we introduce one more. 


Theorem 4.2 (The Resolvent Theorem). Sj k (z ) = —(out/- | (Uo — zl) 1 \irij) or stated 
differently S T (z) = — (Uq — z\)~ l . 


Proof By removing the minor of the a element of U, 


I Un — zl\ + a 


uL J ’ fc) - zl 


<j,k> 


= |Un — zl\ + a IUq — zl\ 


^ k ^ we have that Uq ^ — zl 

where the sub-index indi- 




cates the cofactor of the a element (the \inj) (out k \ element) in U, 
f{z) = |U 0 - zl\ 


U k ) 

a 


So, 


9jk(z) = |U 0 - zl\ <jjk> 


U { j k) - zl 


<j,k> 


= |U 0 - zl\ <3 k> , since this cofactor is not a function of a (that row 

and column is removed). It is a known property of cofactors that if B is the cofactor matrix 
of A (that is; every element of B is the corresponding cofactor of A), then | A| A -1 = B 7 . 
It follows that: 


Sj k (z) = - 


9jk {^) 


|U 0 - Zl\ (U 0 - Zl 


1 kj 


= ~ (U 0 - zl 


] kj 


(16) 


/(A) |U 0 — zl\ 

This applies only to those edges that are “prepared” to be inputs and outputs as 
described at the beginning of this section. That is, a needs to be the only element appearing 
in both its row and column. Replacing an arbitrary element of Uo with a destroys the 
unitarity of the time step operator, U. 


Equation 16 says that when you want Sj k ( A), the scattering coefficient between |inl¬ 


and | out k ), you can find it in the \out k )(irij\ element of (Uq — zl) 1 . In other words: 


Sjk(z) = ~{out k \ (U 0 - zl) 1 \inj) 


(17) 


□ 


The above theorem applies to the single runway case as well. 


4.1 Particulars for Multiple Runways 

Once again, define = Uo+ot\inj){out k \ and Cj k (z) = 


U { f k) - zl 


The zeros of gjk(z) are not fixed by the zeros of f(z), the way t' 


= H z ) ( f ( z ) + agjk(z)). 

rere are in the single 


runway case, and are not necessarily outside of the unit circle. We can see an example of 


this in appendix A.5 
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As in thm. | 2 . 1 [ those eigenvalues with modulus 1 correspond to bound eigenstates, and 
are factored out as b(z). However, in the general case we find that gjk(z) and f(z ) may 
have zeros in common. If | ip) is an eigenstate such that (outkl'f’) = (irij\if) = 0, then it is 
independent of a , however it may not necessarily be a bound state since there are other 
runways that it can be “leaking” out of. If an eigenvalue is independent of a then it must 
be a common factor of /(z) and gjk(z) and if the associated eigenstate has a component 
on any of the | out) states, then its eigenvalue must be less than 1. In fact, in a derivation 
nearly identical to that found in theorem | 2 . 1 | we find that the eigenvalue, g, of a eigenstate 
of U 0 , | if), satisfies \g\ 2 = 1 - YkLi \{out k \^)\ 2 < 1 

These are the zeros of f(z). If (outk\if) = 0, then g is also a zero of gjk(z) and \tf) is 

an eigenstate of UVj. 


Theorem 4.3. If \z\ = 1, then Yk \^jk( z )\ 2 = 1- 

Proof The A-eigenstate for a graph G attached to M runways with a signal coming in 
from the jth runway takes the form 


M 


I*> = £ a ” + 1 I" + 1 -">j+E 

n=0 k= 1 . 


Sjk{ A)^A n \n,n + l) k 


n =0 


M 


+ \G) + \inj) + ^2 XSjk(X)\out k ) 


k= i 


The subscript on the runway states indicates to which runway they correspond. This is 


an eigenstate, so U|\k) = A|T) and it follows that U (|G) + | inj}) = X ^|G) + Yk=i ^^jk 
Being unitary U is an isometry and therefore 
||(j) + \irij)\ = A|G) + A 2 Yk=i Sjk(ty\outk) 

=*■ 1 + <G|G) = |A| 2 (G|0) + |A|‘E", |S, t (A )| 2 

Clearly if |A| = 1, then the statement of the theorem follows immediately. 


It may seem worrisome that A isn’t assumed to be modulus 1 despite being the eigen¬ 
value of a unitary operator, but keep in mind that |\k) isn’t normalizable and U is operating 
on an infinite dimensional Hilbert space. However, on any finite state we can still make 
use of the fact that unitary operations are isometries. 

□ 


4.2 Arbitrary Signals for Multiple Runways 

As before, define the input x[n] (output y[n\) as the amplitude on the state 11,0) (state 
10,1)) at time step n (time step n + 1). The input can be encoded onto the runway in an 
initial state of the form x [ n ]| ra + 1, n). If x\n] = X n for all n, then y[n\ = Sjk(X)x[n]. 


21 






Applying exactly the same proof used in the single input/output case (thm. 2.9) we 
find that 


hjk\n] = 


2ni 


1 


z n L Sjk(z) dz = — -—: 

Zm 


z"-' 9 -4$dz 


fU) 


(18) 


'|a|=l Z7rl J\z\=l 

where hjk[n\ is the response produced by the kth output to an impulse received from 
the jth input. 


Like the single runway case we find that Sjk(z) = Yl^Lo hjk[n]z 
gives us a second proof and a little insight into theorem |4.2| 


(see eq. 12). This 


In general, the impulse response is hjk[n] = (0,1| (U^' fc )) n+1 11,0) (see eq. 11), where 
uO'fc)|i, 0 ) = I irij) and U^| outk) = 10, 1) . The way we have defined the graph (such 
that it includes “in” and “out” states) implies that hjk[0 ] = hjk[ 1] = 0, and therefore 
Sjk(z) = hjk[n]z~ n - In what follows we’ll keep hjk[ 1] in the sum; this changes 

nothing but makes the derivation a little smoother. 


Sjk(z) 

= E”=i h jk [n\z~ n 

= E“,i^“(o.i|(u w )" + 1 iro) 

= Y,n=l z ~ n ( out k\U\ 


n— 1 1 


irij 


= J2n=o z n (out k |Uq | irij ) 

= (out k \ [En=0 Z_n_1U o] I in J 


= (outk | [iEn =0 Q U o)”] |*n. 


= (outk | 
= (out k \ 


) (I - lUo) 
- (Uo - zl) 


-i 


-l 


\irii 


\irid 


So, we can either think of theorem 4.2 as being a result of the nature of the characteristic 
polynomial of U a or as a symptom of the fact that the power series of the negative resolvent, 
— (Uo — ’I) -1 , is a sum of impulse responses multiplied by z~ n which is equal to the 
frequency response, S(z). 
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A Appendix 


The goal of this appendix is to present a series of elucidating examples and to demonstrate 
and verify the theorems presented in the paper. 

The bolo graph example is a single runway graph that’s simple enough for the theorems 
described in this paper to be verified directly through brute force calculation. 

The star graph and complete graph examples reiterate some of the basic theorems from 
section 2 as well as demonstrate how the structure of a graph can be discerned using signals. 

The simple valve and square junction examples are multiple-runway examples. The 
simple valve has two runways and we see how we can either attain total reflection or 
transmission by changing the graph or the signal. The square junction has four runways 
and is included to show how easily the techniques can be generalized. 

The pruned tree example shows how, in practice, the scattering coefficients of a sub¬ 
graph can be used to replace it in a larger graph. 

A.l Example: The Bolo Graph 

The bolo graph is among the smaller non-trivial graphs. It has one bound eigenstate and 
is named for its resemblance to a bolo tie. 


l 



Figure 3: The bolo graph attached to an infinite runway. Here | in) = 10, B) and | out) = 
\B, 0). That is, XJ| 1, 0) = 10, B) and U|F>, 0) = 10,1} whereas U a | B, 0) = a|0, B). 


The states on the Bolo graph are 
IV’i) = \0,B) = I in) 

IV>2 ) = \A,B) 

\ip3) = \B,B) 

1^4 ) = \B,A) 
iV’s) = \B,0) = | out) 

and U a is defined to act on these states as 
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u„ = 


/ 0 0 0 0 a \ 

0 0 0 -1 0 

2 2 _1 Q Q 

a 3 3 U U 

5 1 0 0 

i oo/ 


v A 


2 

3 3 


The characteristic polynomial, written using the form described in theorem 2.1 


is 


|U Q - zl\ = - (z + 1) 

b(z) 


z — — z T — ) T ol [ —z — — z -f- 1 


/(*) 


ffO) 


In this form it’s clear that f[z) and g{z) are reciprocal polynomials. We quickly find 
that rj = | l+h/2, 1- 3 ^ |; s = 2, and go = 1- 

This has five roots and when |ck| = 1 all five have modulus 1, as they should since U a 
is unitary when |a| = 1. r/| = 


lzbi\/2 

3 


= so the roots of f(z) all fall inside the unit 

circle (thm. 2.1). 

Graphing the magnitude of the relevant part of the polynomial (ignoring the “z + 1”) 


demonstrates the relationship between a and the zeros as described by theorems 2.5, 2.6| 


and 2.7 (see fig. 4). 


According to theorem 2.2 the effective reflection coefficient off of vertex 0 is 

z 2 — 2z + 3 


S(z) = - 


9{z) = _ 

f(z) z 2 (3z 2 — 2z + 1) 


(19) 


We can find S(X) by brute force one eigenvalue at a time by calculating the appropriate 
eigenstate. For A = i, we find that the eigenstate is: 


OO OO 1 I ' 1 ' 


3=0 


3=0 


(20) 

Comparing carefully with the definition of S( A) described in section 2, we find that 
S(i) = i. This lines up exactly with what we should expect, since (from eq. 19) S(i) = 
i 2 —'2i+3 


i 2 (3i 2 -2i+l) 


= l. 


Often we would prefer not to have an extra edge on the runway. In this case that edge 
is between B and 0. In order to find the effective reflection coefficient of vertex B instead 
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Figure 4: \f(z) + ag{z)\ for the bolo graph where a = e lx and z = e iy . The zeros are 
clearly visible as troughs. Notice also that for a given zero there is one corresponding value 
of a. 


of vertex 0 we need to delay the signal by two time steps. This is accomplished easily 
through a multiplication by z 2 . The effective reflection coefficient of vertex B is 


R(z) 


2 S(f) __ 


z 2 — 2z + 3 
3z 2 -2z + l 


( 21 ) 


Although the bolo graph only has one input and output state, we can still find the 
scattering/reflection coefficient using the resolvent, as described in theorem 4.2. Using the 
same basis states as before, we find that 


/ 


U n = 


0 

0 

2 

3. 


V -i 


0 

-1 

0 

0 

0 


0 \ 
0 
0 
0 
0 / 


Here the \in)(out\ = |V ; i)(V ; 5| element has been marked in red. We find that the negative 
resolvent is 
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A 


A 2 

- > - * 

• - •-•- 

2 * - 1 * - 0 

A 1 S(A) S(A) 



A 



Figure 5: S( A) is defined normally, as the effective reflection coefficient of vertex 0. R( A) 
is here defined as the effective reflection coefficient of vertex B. From the diagram it’s 
straightforward to see that R(z) = z 2 S(z). 


- (U 0 - ^I )" 1 = 


V 


1 

z 

_ 2 

z(l-2z+3z 2 ) 
2(—1+z) 
z(l—2z+3z 2 ) 

2 

l-2z+3z 2 

3—2z+z 2 
z 2 (1-2z+3z 2 ) 


z(l+3z) 

l-z+z 2 +3z 3 

2z 

l-z+z 2 +3z 3 
1 —z 

1—z+z 2 +3z 3 
2 

1—2z+3z 2 


l-z+z 2 +3z 3 

— l+3z 2 

1—z+z 2 +3z 3 

2z 

1—z+z 2 +3z 3 
2 (—1+z) 
z(l-2z+3z 2 ) 


l+3z 

1—z+z 2 +3z 3 

_ 2 

1—z+z 2 +3z 3 

z(l+3z) 

1—z+z 2 +3z 3 


o \ 
0 
0 
0 


The scattering coefficient is the \out)(in\ = |V^s)(V 2 ! 


z(l—2z~h3z 2 ) z j 

element and is again marked in 


red. As expected (tlrnr. 4.2) this is precisely the result derived earlier in this example. 


Finally, we look at the impulse response of the bolo graph. The first few terms of the 
impulse response can be found by direct calculation: 


n 

U n |l,0) 

x\n] 

y M 

0 

|i,o> 

1 

0 

1 

|0 ,B) 

0 

0 

2 

§| B,B) + || B,A) - l\B,0) 

0 

i 

3 

3 

-11 A,B) - || B,B) + l\B,A) + - ||0,1) 

0 

4 

9 

4 

-II A,B) - %\B,B) + §f\B,A) - $\B,0) + | 0,1) - ||1,2) 

0 

16 

27 

5 

~§f\ A >B) - B,B) - B,A) - ff|5,0> - |0,1) + ||1,2) - ||2,3) 

0 

44 

81 
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Keep in mind here that y[n\ = S(A)a;[n] when x[n\ = \ n . As a result of this definition 
y[n\ “anticipates” when you write out a list like this: during the nth time step cc[n] is the 
coefficient of 11,0) and y[n — 1] is the coefficient of 10,1). 


Rather than an infinite brute force calculation, we can apply theorem 2.10 Plugging 
in r/i = 1+ g v/ ^ , r /2 = 1 ~g V ^ , s = 2, and go = 1 we find: 

On = - 


Vl+*V2 


Oi — — 


Oo — — 


1- 

I+ia/2 

3 

( l+i\/2l 


3 ) 

i- 

l-iV2 

3 

(l 

-iV2\ 



So for n > 2, 

h[n\ = -3% - 2] - i3V2 (^4^)" + * 3 >/2 


A.2 Example: A Star Graph with Differently Marked Edges 

The only information we can ever hope to gain from a signal are the zeros of f(z) and the 
value of go- Every other quantity described so far can be determined from these. 

However, a graph with d edge states can have any of a huge number of arrangements, 
and each vertex in that graph can have any appropriate unitary operator. It is completely 
unreasonable to hope for a way of determining the structure of a completely arbitrary 
graph by studying its responses to signals. That is to say, using the techniques described 
in this paper we cannot hope to distinguish between isospectral graphs, and therefore 
cannot expect to uniquely determine the structure of an unknown graph strictly from its 
response to signals. 

However, if we restrict the graph to being one of a restricted family of graphs, then we 
can expect to get some information about the structure. First we investigate a star graph. 


Assume that we know that there are N + 1 edges radiating from a central vertex c 
which is diffusive (r = — 1 + for every edge and t = -^Aj- between every pair of edges). 

The terminating vertex of some of M of these flip the phase (multiply by -1) when 
they reflect, N — M of them leave the phase unchanged (multiply by 1), and the last edge 
connects to the runway. The question is: can we determine the number of each type of 
edge? 

Define the states 
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IVh) = |0,c) = I in) 

|^ 2 ) = |c, 0) = | out) 

\rk) = jsX?=i\ W) 
1^5) = Ef=! M M. 


V 


s/N-M ^J= 

\^) = 7WTT^U M M 


Here the aj vertices reflect with —1 and the bj vertices reflect with 1. The time step 


operator can be written 
/ 0 - 1 + 2 


U„ = 


N+l 


a 

0 

0 

0 

V 0 


0 

o \Z~M 
A N+l 


o ^N-M 
Z N+l 


0 

0 

0 

-1 

0 

0 


) Vm 

J N+l 
0 


-1 + 


2 M 
N+l 


0 

f y/M(N-M) 
‘ N+l 

0 


0 

0 

0 

0 

0 

1 


r,+N-M 
Z N+l 
0 

, y/M(N-M) 
J N+l 
0 


-1 + 


2(N—M) 

N+l 

0 


The characteristic polynomial is 


C(z, a ) = z 6 + ( a — 


2 a 


+ 


AM 


2 N 


N+l N+l N+l 


z 4 + -1 + 


AM ot 


+ 


2N 


2 Not 


N + l N + l N + l 


z A +a 


= z 


z* + 


AM - 2N 
N + l 


z z + 


N -l 
l+N 


T OL 


N - 1 
N + l 


z* + 


AM - 2N 
N + l 


z 2 + l 


m 


s(*0 


The effective reflection coefficient of vertex c is 

Rlz] = z 2a,s = (JV-l) Z 4 +(4M-2Af)^ + (JV + l) 

y ' K > (N + l)z‘ + (4M - 2JV) 2 2 + (jV - 1) 


Where the extra z 2 is a delay necessary when the 2 extra states, |0,c) and |c, 0), are 
included in the runway. It is not immediately obvious how this function behaves. However, 
we can hope to understand it better by looking at the zeros of f(z) and applying theorem 


3.1 The four non-zero roots of f(z) are 


T] = + 



2M + 1 y/1 - AM(N - M) 

N + l ± N + l 


( 22 ) 
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where ± and ± are independent. 


Define L = then 

V = ± ^(1 " 2L) (1 - £) ± 2iy/L(l -L)(l-±)+0{jk) 

= ±(1 - ^)y/(l-2L)±2iy/L(l-L) + O (^) 

= M 1 ~ e±i2T + ° iw) 

= ^( 1 ~ m) e±lT + O iw) 

where sin(2r) = 2y / L(l — L). 

These four roots are negatives and/or complex conjugates of each other. The problem 
of finding L = ^ can now be reduced to finding the phase of the root of f(z) in the first 
quadrant, which is made difficult by the fact that rj « (1 — 7^)e* r lies just inside of the 
unit circle. Away from these zeros the frequency response is fairly unrevealing. 



Figure 6: Re [i? (e* 0 )] for N = 100, M = 40. 

With the exception of the thin spikes that occur around the four zeros of f(z), this 
approximately equal to — 1, which is the frequency response of a reflection by — 1 at vertex 
c (which is approximately what the graph is). 

This means that in order to detect one of the r/'s we need to detect an interval in which 
S(e le ) 7^—1, and that difference is only detectable within () (^) of the correct value. 

So, if a pulse is sent and it is found that S(e ld ) / —1, then 2y / i?(l — R ) = sin(2$). 
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A.3 Example: Complete Graph of Unknown Size 

Here we consider a complete graph with N + 1 vertices, labeled Aq,Ai, ..., Ajv- All of 
these are connected to each other, and A$ is additionally connected to vertex 0. We define 


U| A jt A k ) 



I Afc,A ? -) + — ^ \A k ,Ai) 
i=o 


and 


ui A hA)) - (-1 + jvq^;) l^’^ + jv+T 2 + ivTT 

v 7 i=i,i & 


1^0,0/ 


Define the basis states as: 

IVh) = |0, A 0 ) 

\i’2) = ^EU\ A o,A j ) 

^ = y/N(N-l) =1 'Ek = l,kft \ A Jl A k) 

l^) = ^Ef=i|A i ,A 0 ) 

IV's) = |A 0 ,0) 


We find that U a is: 

/ 0 


U„ = 


2 VN 

JV+1 

0 

0 


V -1 + 


0 

0 

2 VF 4 

N 

-1 + A 

1 ^ N 


0 

0 


1 - A 

1 TV 

2WV-T 

N 


1 - 


JV+1 


JV+1 


0 

0 

2 yW 
JV+1 


a \ 
0 
0 
0 
0 


It follows that the characteristic polynomial of is 

r ,, \_ 2 (N 2 +N)z 3 +(2+N-N 2 )z 2 +(2-3N+N 2 )z+(N-N 2 ) (N-N 2 )z 3 +(2-3N+N 2 )z 2 +{2+N-N 2 )z+(N 2 +N) 

<-y{Z,a)—Z _ JV(JV+1) a N{N+1) 

and therefore (by thm. 2.2) the effective reflection coefficient of Aq is 


(N - N 2 )z 3 + (2 - 3 N + N 2 )z 2 + (2 + N - N 2 )z + ( N 2 + N) 
(N 2 + N)z 3 + (2 + N - N 2 )z 2 + (2 - 3A + N 2 )z + (N - N 2 ) 


Again, S(z) is the coefficient for vertex 0 and R(z) is the coefficient for vertex Aq. 
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Like in the last example this equation couldn’t be less clear, but like the last example 


we can gain some insight by using thm. 3.1 By writing f{z) as 


f(z) = z 2 (z 3 - (1 - #) - 2 + (l + jf ~ atti) z ~(*-n+ T y 

= z 2 (z 3 - z 2 + z - 1) + Jr (2z 2 -4z + 2) + j^r(6z-2)+ 0 ( N~ 3 ) 

= z 2 (z 2 + l)(z-l) + ^(z-lf + ^ (; 3z - 1 ) + 0 (IV- 3 ) 
we can quickly find that to lowest order the roots of f(z) are: 

1 + * . 1 — 1 2 


0 , 0,1 — 


IV 


1 1 — 1 

-•‘'-IT’ 1 - 


jv 2 


The last of these is the most difficult to detect and has a constant phase, which means 
that it is a poor choice for a probing signal. The other two roots do change phase (a little) 
and are easier to detect. We find that, again to first order, arg (l — ^) = f + and 
\i — 1+i I = i_L 

U JV I 1 AT 


N ' 



Figure 7: Re [i? (e 20 )] for N = 10. The vertical lines indicate 6 = ±|. 

Once again, this effective reflection coefficient of Aq is approximately —1 for most values 
of 6, which is to be expected since (Ho, 0|XJ|0, Aq) = —1 + qqq- ~ —1. However, if we can 
find the value of 6 such that R (e* e ) = 1, then we can determine the number of vertices in 
the complete graph since R (e* e ) = l=># = ±(f + -^)+ 0{N~ 2 ). 
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A.4 Example: A Simple Valve 

Here we have a degree 3 vertex, D , connected to three other vertices: A, B , and C. D is 
a standard diffusive vertex, so t = | and r = — A and B reflect with 0 and C reflects 
with some constant c. When attached to a runway vertices A and/or B are vertex zero. 



Figure 8: The simple valve. When attached to a runway A and/or B will be vertex 0. 


Define the basis states as 


IVh) 

= 1 A,D) 

= \in\) 

1^2) 

= 1 B,D) 

= | m 2 ) 

l^3> 

= 1 C,D) 


1^4) 

= 1 D,C) 


1^5) 

= 1 d,b) 

= \out 2 ) 

IV’e) 

= 1 D,A) 

= \outi) 

The time step 

operator, 


/ 0 0 0 0 O 12 On \ 

0 0 0 0 O22 O21 

0 0 0 c 0 0 

u "= I § -1 0 0 0 

I -s 1 0 0 0 

— - - -ODD 

\ 3 3 3 U U 

In the above matrix the locations corresponding to the terms that scatter \irij) into 
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\outk) are marked in red and labeled 0 jf.. Each of these entries are in fact zero, merely 
labeled. Were we using thm |2.2| (instead of thm 4.2), then 0 jk is where a would be inserted 
in order to calculate Sjk{z). 

A quick (computer aided) calculation reveals that 

( 


-(Uo-^I )- 1 = 


V 


1 

z 

0 

0 

0 

0 

0 

0 

1 

z 

0 

0 

0 

0 

2c 

2c 

3 z 

3c 

0 

0 

z(3z 2 +c) 

z(3z 2 +c) 

3 z 2 +c 

3z 2 +c 

2 

2 

-1 

3z 

0 

0 

3 z 2 +c 

3z 2 +c 

3z 2 +c 

3z 2 +c 

2 (z 2 +c) 

-z 2 +c 

2 

2c 

1 

0 

z 2 (3z 2 +c) 

z 2 (3z 2 +c) 

3z 2 +c 

z(3z 2 +c) 

z 

-z 2 +c 

2 (z 2 +c) 

2 

2c 

0 

1 

z 2 (3z 2 +c) 

z 2 (3z 2 +c) 

3z 2 +c 

z(3z 2 +c) 

z 




where again, the relevant terms are marked in red. 

A subtle point to keep in mind is that these describe 11,0) scattering into |0,1), which 
are not included in G. To find the effective reflection coefficient of vertex D we simply 
need to get rid of the delay caused by the two extra edges, and this is done by multiplying 
by z 2 (see fig. [HJ) . 

We see immediately that the effective reflection and transmission coefficients (from 


either direction) of vertex D are r +> = w and *(*) = respectively. 


When z 2 = c we see that D becomes “transparent’' 

2(c+c) 


since r(± N /c) = - 


= 0 and 


C+C 


*(±VS) = = 1- 

When z 2 = —c we see that D becomes completely reflective, since r(±iy / c) = + c+c — 
-1 and t(±iy/c) = = 0. 

This also provides a tool for determining the value of c: if a particular momentum state 
with eigenvalue A is entirely transmitted, then c = A 2 . In fact, this method was used in [8] 
to distinguish between the two results of a boolean function. When the effective reflection 
coefficient of a tree graph constructed to execute a particular calculation and attached to 
vertex C was —1 (1) the ii-eigenstate was entirely transmitted (reflected). 


A.5 Example: Square Junction 

In this example we look at a graph with four inputs and outputs and discuss what the 
scattering coefficients mean in terms of eigenstates and signals. Even with this (16 state) 
graph we begin to see the use for replacing graphs with a single vertex (with non-trivial 
scattering coefficients). 

Define the basis states as 
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Figure 9: A four-way junction. 


IV’i) = k A) = I ini) 
k2> = | b,B) = | in 2 ) 
IVk = I c,C) = | m 3 ) 

| ^ 4 ) = k-0) = Kn4) 

m = \A,B) 

\A) = \ B,C) 

Vh) = \c,d) 

ks) = 1-0, A) 
k9> = l A,D) 
kio> = | D,C) 
kn) = l C,B) 
ki 2 ) = l-B, A) 
kis) = 1-0, d) = |out 4 ) 
ki4> = | C,c) = | out 3 ) 
kis) = 1-0,6) = | out 2 ) 
Ike) = I A, a) = \outi) 
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The time step operator, Uq, is 


0 
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0 

0 
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021 
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034 
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044 
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0 

0 
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0 

0 

2 

3 

0 

0 

1 

3 

0 

0 

0 

2 

3 

0 

0 

0 

0 

0 
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0 

2 
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0 

0 

1 
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3 
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0 
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0 

0 

0 

1 

3 

0 

0 

2 

3 

0 

0 

0 

0 

0 

2 

3 

0 

0 

0 

0 

0 

1 

3 

0 

0 

0 

0 

0 

0 

2 

3 

0 

0 

0 

2 

3 

0 

0 

0 

0 


The gridlines here are included to make the matrix more intelligible. The zeros in the 
upper-right are labeled according to which input and output they correspond. 

To find all of the scattering coefficients we use thm. 4.2 and find that the 16 lower left 
most entries of the negative resolvent are 


-(Uq-^I )- 1 



\ 

4 

16 

4 

9z 4 +10z 2 —3 


z(9z 2 -1) 

27z 4 +6z 2 —1 

Id 

bi 

to 

1 

z 2 (27z 4 +6z 2 —1) 


16 

4 

9z 4 +10z 2 —3 

4 


27z 4 +6z 2 -l 

^D 

10 

z 2 (27z 4 +6z 2 —1) 

z(9z 2 —1) 


4 

9z 4 +10z 2 —3 

4 

16 


ts* 

Id 

to 

1 

z 2 (27z 4 +6z 2 —1) 

z(9z 2 —1) 

27z 4 +6z 2 —1 


9z 4 +10z 2 —3 

4 

16 

4 


z 2 (27z 4 +6z 2 —1) 

81 

'CD 

10 

1 

27z 4 +6z 2 —1 

z(9z 2 —1) 

/ 


/ 

(23) 


We can gain a little more insight by instead applying thm. 2.2 to find S'n (z). First we 


look at the characteristic polynomial obtained by replacing On with a: 
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Cii{z) = z 8 (z 8 - lz 6 - |fz 4 - ^z 2 + ±) + az 6 (±z 8 + ±z 6 - jjfz 4 - ±z 2 + 

= || (3 z 2 + l) ( z 2 — l) [z 2 (27z 4 + 6 z 2 — 1) + a (9 z 4 + 10 z 2 — 3)] 

In this form we see why the scattering coefficients takes such a simple form: in each 
case f(z) and gjk{z) have several common factors. 

The terms that can be factored out of both f(z) and gn(z) correspond to eigenstates 
that are uninvolved with signals scattering between \in\) and |crati). For example, there 
are two bound states with eigenvalues ±1 that aren’t involved in any signals (being bound, 
they have no overlap with any of the input and output states). 

Since 5] 1 (z) = — , the zeros of gu{z) indicate which signals from \in\) will return 

nothing to \out\). These zeros are A = |±|y / 2\/l3 — 5, V / 2\7l3”-|-~5 j. This means that 

a signal of the form x\n] = (j^\/2y/ 13 — coming in from any of the inputs will have 
no reflection; it will scatter entirely out of all of the other outputs. It’s worth pointing out 
that, in this example, gjk(z) is trivial when j ^ k. As a result every signal will scatter out 
of all of the other outputs at least a little (since gjk(z) has no zeros). For example, there 
is no signal that can be sent through \iri 2 ) that will not produce an output through \ 0 ut 3 }. 

Notice also that the real roots are inside of the unit circle and the imaginary roots are 
outside. This is a beautiful example of the zeros of gjk(z) neither being determined by the 
roots of f(z) nor being strictly outside of the unit circle (as they are in the single runway 
case). 


By using these 3 equations (eq. 23) for the square junction we can easily insert this 


graph into a larger one without adding an additional 16 states, which is important if you 
want to keep your computational overhead low. 


A.6 Example: The Pruned Tree 

One of the central claims of this paper has been that a subgraph can be replaced with a 
single vertex with a scattering coefficient we can calculate. In this example we demonstrate 
this claim. 

We calculate the scattering coefficient of a two-level tree, G, by two methods. First, we 
calculate it directly using the theorems established in this paper. Second, by replacing the 
second level of the tree with a scattering coefficient, r(z), reducing it to a one level tree. 

These two methods yield identical results and I argue are in fact equivalent. That 
is, replacing a subgraph with it’s scattering coefficient is a useful way of simplifying and 
understanding the calculation of the overall scattering coefficient. 

We define the basis states as follows 


39 








I 

I 


Figure 10: “The Pruned Tree”. We can calculate the scattering coefficient of the graph as 
a whole by replacing G with G', and G" with the reflection coefficient r(z). 


IVh) = 

IA0) 

IV>2> = 

|0,^> 

1^3} = 

1 A,B) 

1^4} = 

\B,A) 

1^5) = 

1 A,C) 

\A) = 

\C,A) 

IV't) = 

1 D,C) 

IV’s) = 

1 C,D) 

ivto) = 

1 E,C) 

IVho) = 

- \C,E) 


and order them such that the full graph is G = {|V’i} ) • ■ •, IV'io)}s the simplified graph 
is G' = (IV’i), • • •, |-06 )}5 and the subgraph to be “pruned” is G" = (l^),..., IV'io)}- 

First we find the scattering coefficient of the entire graph. 
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The full time-step operator is 
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0 

0 

0 

0 

-1 
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0 

0 
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2 

3 

0 

2 

3 

0 

1 

3 

0 


The characteristic polynomial is 

C{z) = z 2 (z s + §£ 6 + %z 4 z 3 + l)+a (lz 8 - |^ 6 + %z 4 + \z 2 + 1 ) 

and according to theorem 2.2 the scattering coefficient of vertex 0 is 


S(z) = 


3z 8 - 2 z 6 + 4z 4 + 2z 2 + 9 
z 2 (9z 8 + 2 z 6 + 4z 4 - 2z 3 + 3) 


This is the reflection coefficient for the entire graph. We want the results derived below, 
where we’ve “pruned” the graph, to match this. 


The subgraph uses the basis states |^) through IV’io) » and has the time-step operator 


Q 

o 

0 0 0 0 \ 

1 

o 

2 0 - 0 

3 u 3 u 

0 0 

0 10 0 

2 0 

3 u 

3 u 3 u 

0 0 

0 0 0 -1 

2 0 

3 u 

1 0 -3 0 ) 


The characteristic polynomial of this is 
C"(z) = z 2 {z 4 + \)+a(\z 4 + 1) 
and again by theorem |2.2[ 


S"{z) = 


z 4 + 3 
z 2 (3^ 4 + 1) 
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To find the correct scattering coefficient, the reflection coefficient of vertex C, is 


r(z) = z 2 S"(z) = 


z 4 + 3 
3z 4 + 1 


Again, this z 2 is to adjust for the A, C edge (see fig. [5| 


The “pruned tree graph”, G ', uses the basis states through \^%). We replace the 
entire subgraph, G", with a reflection coefficient r(z) such that U'l A, C) = r(z)\C, A). We 
can use this when calculating eigenstates, and thus when calculating scattering coefficients. 
The pruned tree graph is G' = ..., \ipe)} and its time-step operator is 


0 4 0 i 

0 1 

o 

o 

o 

Q 

0 0 

0 - 0 —- 

u 3 u 3 

0 2 

u 3 

0 0 10 

0 0 

0 I o I 

0 -3 

\ 0 0 0 0 

r(z) 0 


The characteristic polynomial is 
C'(z) = z 2 


z 4 + i±rM z 2 _ 


Kf) 

3 


+ a ( T^z 4 — 


l+r(z) „ 2 

Q .<✓ 


-r(^)) 


The reason this should make sense is that the effective reflection coefficient is the 
response that graph has as part of an eigenstate. Indeed, if A is an eigenvalue of the full 
graph, then it should satisfy C'{ A) = 0. From theorem 
coefficient of the graph is 


2.2 


we know that the scattering 


, z 4 - (1 + r(z))z 2 - 3r(z) _ z 4 - z 2 - [z 2 + 3) r(z) 

z 2 (3z 4 + (1 + r(z))z 2 — r(z)) z 2 (3z 4 + z 2 + {z 2 — l)r(z)) 

If we correctly represented the subgraph, then when we insert r(z) into S'(z) we should 
recover the full S(z). 

S'(z) 

_ z 4 -z 2 -(z 2 +3 )r(z) 

z 2 (3z 4 +z 2 +{z 2 -l)r(z)) 

= z 4 -z 2 +(z 2 + 3 ) ^ 

Z 2 (3Z 4 +Z 2 -(z 2 -1)|^L) 

_ (z 4 -z 2 )(3z 4 +l) + (z: 2 +3)(z 4 +3) 

~ ~ z 2 ((3z 4 +z 2 )(3z 4 +l)-(z 2 -l)(z 4 +3)) 

_ _ 3z 8 -2z 6 +4z 4 +2z 2 +9 

z 2 (9z 8 +2z e +4z 4 -2z 3 +3) 

= S(z) 
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So, the scattering coefficients obtained from the pruned tree, G', and the fuli tree, G , 
are the same. 

The four poles of r(z) are of no great concern. As S'(z) = S(z ) for all z for which r(z) 
is defined, these are removable singularities. 

So, the scattering coefficient of a graph can be found by replacing a subgraph with its 
own scattering coefficient. This is useful if, for example, you wanted to do a calculation 
involving a graph with a lot of identical subgraphs. 
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